Binding and structure of tetramers in the scaling limit 
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The momentum-space structure of the Faddeev-Yakubovsky (FY) components of weakly-bound 
tetramers is investigated at the unitary limit using a renormalized zero-range two-body interaction. 
The results, obtained by considering a given trimer level with binding energy B3, provide further 
support to a universal scaling function relating the binding energies of two successive tetramer 
states. The correlated scaling between the tetramer energies comes from the sensitivity of the 
four-boson system to a short-range four-body scale. Each excited N— th tetramer energy B^ 
moves as the short-range four-body scale changes, while the trimer properties are kept fixed, with 
the next excited tetramer B^ N+1 ^ emerging from the atom-trimer threshold for a universal ratio 
B4 /B3 = B^ NS> I B^ N+1 ^ ~ 4.6, which does not depend on N. We show that both channels of the 
FY decomposition [atom-trimer (K— type) and dimer-dimer (H— type)] present high momentum 
tails, which reflect the short-range four-body scale. We also found that the //—channel is favored 
over K— channel at low momentum when the four-body momentum scale largely overcomes the 
three-body one. 
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I. INTRODUCTION 

A limit cycle in physics [1] refers to a model indepen- 
dent way to look for a manifestation of a hidden scale, 
one in which a correlation between physical quantities is 
geometrically rescaled and replicates itself. Limit cycles 
are shown to exist even in simple quantum three-boson 
systems [2H3], manifested for large two-body scattering 
lengths, as well as for short-ranged interactions. They 
follow the 1970 Efimov key prediction [5] of geometri- 
cally separated weakly-bound three-boson states near the 
scattering threshold, confirmed by recent experiments in 
cold-atom traps [SHE]- For a more recent exposition on 
the experimental status, see Ref. [§]. The addition of 
one more particle to the quantum three-body system has 
long challenged this picture [TDHT5] . Recent theoretical 
works [T4Tf2"T] motivated by cold-atom experiments [%2j - 
[53] have revived this issue. The actual relevance in identi- 
fying simple universal relations in few-body binding laws, 
with the corresponding experimental possibilities, have 
been discussed by Modugno 25i in a recent perspective 
article on this matter. 

In Ref. [25], we reported a different dimension of 
this problem by establishing a new universal correlation 
among the binding energies of two successive tetramer 
states between two Efimov trimers, through precise nu- 
merical calculations within a zero-range interaction at 
the unitary limit (zero two-body binding). It was 
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shown that such tetramer states are related to an un- 
suspected new limit cycle not constrained by Efimov 
physics or by three-body properties. Tetramers hitting 
the atom-trimer threshold, leading to resonant recombi- 
nation losses in this channel, are signatures of a new limit 
cycle. Furthermore, other four-boson observables close to 
the unitary limit can exhibit correlations not constrained 
by the low-energy two- and three-boson properties. Re- 
combination rates, atom-trimer or dimer-dimer scatter- 
ing lengths, can move near the Feshbach resonance (FR) 
independently of the two and three-body properties. The 
appearance of few-body forces in the open channel, due 
to the coupling with the closed channel, can drive the 
four-body physics |15) . Indeed, there is a recent exper- 
imental evidence, reported by Nakajima et al. [57], of a 
three-body force acting in the open channel near the FR 
in a three-component mixture of a 6 Li cold gas. 

An appropriate way to directly probe the universal 
properties of a few-boson system at low energy near the 
unitary limit is to present results in terms of scaling func- 
tions correlating pair of observables, within a renormal- 
ized zero-range two-body interaction. In this way, no 
effect can be claimed to be originated from a particular 
form of the short-range interaction between the particles, 
as evidenced by the three-body universal scaling function 
derived in Ref. [2] ■ As shown in this case, the regulariza- 
tion parameter in the renormalized zero-range approach 
can be directly associated to a three-boson physical in- 
formation, in addition to a two-body observable, such as 
the scattering length. The concept of scaling functions, 
which is used to evidence universal properties of three- 
boson systems close to the unitary limit, is extended to 
four-boson systems in Ref. [26], as well as in the present 
work. Within such purpose, we study the sensitivity 
of tetramer properties to a four-boson scale, by using 
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the Faddeev-Yakubovsky (FY) formalism [251 HH] with 
zero-range two-body interactions. Therefore, indepen- 
dently on the choice of the three-body scale that fixes 
the trimer properties [15] , a four-body regularization pa- 
rameter (four-body scale) is introduced in the FY for- 
malism. The new short-range scale has indeed effects on 
four-boson physics, which appears as universal correla- 
tions between four-boson low-energy s— wave observables 
not constrained by the trimer properties. The four-body 
quantities change as long as the new scale varies. 

Previous conclusions for the nonexistence of a proper 
four-body scale rely on a strong suppression of the short- 
range physics, beyond that already accounted for in the 
three-boson system. Indeed, by estimating the trace of 
the kernel of the four-body equation in momentum space, 
Amado and Greenwood [10] showed that there is no in- 
frared divergence, which led them to conclude against the 
existence of Efimov effect in the case of four or more par- 
ticles. However, the momentum integrals should also im- 
plicitly have an ultraviolet cutoff (the four-body one) to 
regulate them. Note that according to Weinberg's eigen- 
value criterion [3U] the divergence of the trace is a nec- 
essary condition for an infinite number of bound states. 
However, in the kernel of the FY equation, this criterion 
can only by applied when considering the trimer in the 
ground state. More explicitly, in case of ultraviolet di- 
vergence of the four-body kernel, by moving a four-body 
scale in relation to the three-body one, an infinite num- 
ber of tetramer states emerges from the threshold, but 
the number is restricted by the energy interval between 
trimers (when the reference trimer is an excited one), 
such that it can be infinite only in case the given trimer 
is the ground state. 

It is also instructive to recall other pioneering four- 
body calculations that have verified the possibility of 
Efimov-like behavior in four-body systems, such as the 
works of Kroger and Perne [11] and Tjon and Naus [13] . 
Conclusions drawn within the nuclear physics context are 
obviously limited, in view of the strong nuclear repulsion 
of the potential core, such that the possible observation 
of a four-body scale, independent from the three- and 
two-body ones, is suppressed. It explains why the 4 He 
and the triton binding energies are strongly correlated 
with a fixed slope, as verified by Tjon |31j . In our under- 
standing of the four-body problem, it is not enough only 
two parameters, which determines trimer properties, to 
describe the four-boson system. Our results imply that 
the correlation between the tetramer and trimer binding 
energies forms a family of Tjon lines with slopes depend- 
ing on the new four-boson scale. To exemplify that, we 
have shown a four-boson universal correlation among the 
energies of two successive tetramers appearing between 
two consecutive Efimov states [33]. The correlation ex- 
hibits a dependence on a new scale not fixed only by the 
trimer properties. In order to directly address this novel 
universal behavior, we have performed a number of calcu- 
lations of tetramer properties within a zero-range model 
to show how the dependence on the new short-range scale 



is evidenced through their structure in momentum space. 

The tetramer energies can be presented in a scaling 
plot where the behavior of a four-body scale can be eas- 
ily verified in comparison with other relevant few-body 
scales [26 . The existence of a short-range four-body 
scale, which is independent of the two- and three-body 
ones, is expressed in terms of scaling functions or cor- 
relation between observables. Once the short-range pa- 
rameters are eliminated in favor of physical quantities, 
the renormalized results appear as correlations between 
the four-body observables with fixed dimer and trimer 
properties. 

In the present work we provide details for the approach 
that lead us to report [23] a universal correlation be- 
tween the energies of two successive tetramers. As it will 
be thoroughly presented in the next sections, we estab- 
lish by our numerical investigations that a new short- 
range parameter is necessary to fix the properties of the 
four-boson system in the unitary limit. Our conclusions 
are supported by precise numerical solutions of the four- 
boson FY integral equations in momentum space for a 
zero-range two-body interaction [15 , as well as by an 
analysis of the four-body wave function. 

The paper is organized as follows: in Sect, [n] we 
present the FY formalism, within a renormalized zero- 
range interaction. The relevant scales are introduced in 
the formalism through a subtractive renormalization ap- 
proach. In Sect. |III[ the scaling functions for the binding 
energies of trimers and tetramers are discussed, follow- 
ing Refs. [2] and [25]. Numerical results close to the uni- 



tary limit are also reported in this section. In Sect. IV 



we present our main results for the structure of univer- 
sal tetramers. We include in this section an extended 
analysis of the low and high momentum structures of 
the atom-trimer [(A+T) or K— type] and dimer-dimer 
[(D+D) or H— type] components of the FY decomposi- 
tion, to pin down where the manifestation of the short- 
range four-body scale is more evident. The behavior of 
the four-body FY components of the wave function is 
presented in this section. Finally, in Sect. [V] we present 
our conclusions with perspectives in relation to possible 
experimental observations. Three appendices are supply- 
ing further details on our approach: In [A] and [B] we have 
details on the four-body formalism and total wave func- 
tion in momentum space. In[C] we discuss stability and 
convergence of our numerical approach, and give some 
details on the Lanczos-type procedure for solution of the 
coupled FY equations. 



II. FADDEEV-YAKUBOVSKY FORMALISM 

The treatment of four-body problems is quite well 
known in the quantum scattering theory, following the 
original Faddeev formulation of the three-body sys- 
tems [28] . later on extended to N— particle scattering 
by Yakubovsky [2J5]. The actually known Faddeev- 
Yakubovsky formalism has been considered by several au- 
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thors in their effort to solve different aspects of four-body 
systems, in nuclear and atomic physics. In addition to the 
works on four-body systems already cited in the introduc- 
tion, we can also mention several other reference works 
concerning details of the bound-state four-body FY for- 
malism, such as Refs. [32H4T] and references therein. Al- 
ternatively, in the solution of four-body systems with 
separable interactions, several authors (exemplified by 
Refs. [42 45 ) have applied the Alt-Grassberger-Sandhas 
formalism 46 , in which the t-matrix components are 
used instead of the wave-functions. 

Here, we describe the FY formalism for bound states 
of four identical bosons, by considering general separa- 
ble two-body interactions. Within our approach, the 
form-factors of the separable potential are conveniently 
replaced by point-like interactions, within a renormalized 
zero-range model where the three- and four-body scales 
are introduced. 



A. Bound-state equations and notation 

In a four-body (4B) system (with particles i, j, k and I) 
there are 18 different coordinate systems, each one asso- 
ciated with a specific two-body partition (see Appendices 
|A| and [b]) . Basically, by considering the arrangements of 
the corresponding two- and three-body sub-systems, two 
different partitions are possible: K— type (A+T), where 
a single particle is bound to a three-body subsystem; and 
H— type (D+D), where two dimers are bound. Clearly, 
as we know from the Faddeev formalism, each three-body 
subsystem has also three atom-dimer partitions. 

We present below the FY formalism with our notation, 
which is used along the next sections. The bound state 
of four particles which interact via pairwise forces Vij is 
given by the Schrodinger equation, 



(i) 



i<j 



where \tpij) = G Vij\ty, G = (E-Hq)^ 1 is the free four- 
body resolvent and Hq stands for the free Hamiltonian. 
The above components satisfy: 



kl^ij 



(2) 



GoUj \ip ik ) + \ipu) + \ip jk ) + \ipji) + \tpki) 



where the two-body transition operator is tij. The FY 
decomposition of \ipij) reads: 



\Ki&) = Go^MVW + l^)), (3) 
\Hij : ki) = Gotij\ipki)i 



y fe ) and \Hij t ki) correspond, respectively, to a 



where \K, t 

A— type and H— type partitions 
nents of the wave function 



The Faddeev compo- 



\Kij,k) 



+ 



ioM) 



(4) 



are reconstructed through FY components. 

Every component contains two K~ type and one 
H— type configurations. Therefore, the total wave func- 
tion | Vt) contains twelve different K— type chains and six 
H— type chains, leading to eighteen independent FY com- 
ponents. By considering identical bosons, the four-body 
wave function |^) has to be totally symmetric. As a 
consequence, all twelve K— type components are identi- 
cal in their functional form with particles permuted. The 
same is true for the six H— type components. Thus, it 
is sufficient to consider only two independent FY compo- 
nents corresponding to the AT— and H— type partitions, 
\K) = \KA k ) and \H) = \H ijM ). The 18 coupled FY 
equations, for identical bosons, shrink to two coupled ho- 
mogeneous equations, 



\ K ij.k) — G tijP 

\Hij,kl) = GotijP 



1 + Pki) \K 



kl I \^ij,kl 



\Hij,ia) 



(•5) 



where Pki is the permutation operator for the pair (kl), 
with P and P defined by 



P = (Py + P ik )P jk and P = P lk P jl . 



(6) 



The symmetry property of \K) under exchange of par- 
ticles i and j, and \H) under separate exchanges of par- 
ticles i,j and k, I guarantee that the full wave function, 



l + P + P k iP + P 



1 + P kl ) \Kil k ) 



\H, 



ij,klj 



l + {l + P)P kl )[l + P)\K i l k ) 



1 + P\\H, 



Hj,kl) 



(7) 



is totally symmetric. In order to get insight on how the 
short-range four-body scale is reflected in the bosonic 
wave function for a zero-range force, we will analyze the 
momentum structure of the K and H components. Fol- 
lowing this strategy we will be able to map where the 
new scale is more relevant to built the tetramer wave 
function. For that purpose, we still simplify the problem 
using s-wave one-term separable potentials. This class 
contains, in particular, the contact interaction. 
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B. One-term separable potential with s —wave 
projection 

In this subsection, to simplify the notation, we use the 
particle-labels 1,2,3 and 4. So, we rewrite Eq.([5]) as: 



\K) = G t 12 P 
\H) = G Q t 12 P 



(l + P 34 )\K) + \H) 
1 + P 34 )\K) + \H) 



(8) 



In order to solve the coupled Eqs. {[sj) in momentum 
space one should project these equations into standard 
sets of Jacobi momenta, corresponding to both, if— type 
(| Ui u 2 u 3 )) and H— type (| vi v 2 v 3 )) partitions, as rep- 
resented in Fig. [TJ 





FIG. 1: (Color on-line) Definition of the four-body Jacobi mo- 
menta corresponding to the if— and H— type fragmentations. 

The standard Jacobi momenta for the 4B system can 
be defined in terms of the single particle momentum vari- 
ables 1^=1^,3,4, in two possible configurations (If— and 
H— type), as shown in Fig. [lj 



Ul 

u 3 

Vl = 

v 2 = 



(1/2) (k 
(1/3) 
(1/4) 



2k 3 
3k 4 



ki 
ki 



(l^Xki-ka), 
(1/2) (k 1+ k 2 ) 

V 3 = (l/2)(k3-k 4 ), 



-k 2 )], 

-k 2 + k 3 )] , 
(k 3 + k 4 )] , 



(9) 
(10) 



where we assume identical particles, with mass m =1. 
Since we are interested in the s-wave channel contribu- 
tion, we introduce the partial-wave representation of the 
four-body projection operators corresponding to each Ja- 



cobi momenta set, as 

\u) = \uiu 2 u 3 ) 



\Vi v 2 v 3 ). 



(11) 



Within the s-wave projection, we consider the following 
completeness relation for both basis sets: 



D 3 U \U)(U\ = 1, 



(12) 



where U indicates each one of u and v sets, with D 3 U = 
Ui dUi U 2 dU 2 U 3 dll 3 . Clearly, the projection operators 
| it) and \v) are adequate to expand |if) and \H) compo- 
nents, respectively. Consequently, the projection of the 
coupled equations (|S|) is given by 



(u\K) = (u\G tP(l + P 34 )\K) 
(v\H) = (v\G tP(l + P 3i )\K) 



(u\G tP\H), 
(v\G tP\H). (13) 



In the following, the units are such that h = 1 and m = 1. 

By considering a one-term separable two-body po- 
tential operator, V = A|x)(x|j the s— wave two-body 
t— matrix elements can be written, in the if— and 
H— representations, by 

WW) = ^ 6 % U2) m ~ U3) xiUiMUl) r(ea), 



(14) 

where the corresponding if— type and H— type Jacobi 
sets, {ui,u 2 ,u 3 } and {v±, v 2 , v 3 }, are represented by 
{Ui,U 2 ,U 3 }. In the above, t(ejj) is the reduced scat- 
tering amplitude, given by 



r(£u) 



d 3 P 



x 2 (p) 



u — p A + id 



(15) 



where, for the if— type and H— type configurations, re- 
spectively, we have, 



f 



E - 



2u 3 , „ 

— - and E V =E 



4(16) 



Note that, £ M and £ v have different expressions in the if— 
and H— representations, but are identical when written 
in terms of the single particle momentum coordinates k;. 

The final expressions for the FY components, derived 
in the Appendix |A| through Eqs. ( A1||A15 1, are: 



iff Hl(u' 2 ,U2,x),U2,U3 ) + 



if (ui, 112,113) = 47r Go 112,2/3) xi u i) T (£u) J du'211'2 J dxx\^li(u2,u' 2 ,x) 
x I if ^IIi [u 2 , «2, x) , n 2 (w' 2 , «3, x') , II3 (u 2 , «3, a:')^ + HUIi (u 2 , M2, a:) , II4 (u 2 ,u 3 ,x'),U 5 [u 2 ,u 3 ,x') 



H(vi,v 2 ,v 3 ) = 4:TvG (v 1 ,v 2 ,v 3 ) x(vi)t(£ v ) / dv' 3 v'£ x(4) 



dxK I v 3 ,U 6 (v 2 ,v' 3 ,x) ,II 7 (v 2 ,v 3 ,x) ) + H ( v 3 , v 2 , v' 3 



dx 



(17) 
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where the functions IIj— 1,2. ..7 are defined by the 



Eqs. (A5l, (A10), (All), (A13I, (A14), (A16) and (A17l 



Eqs. (17) can be rewritten within new definitions for the 



FY components, such as 



K(ui,u 2 ,u 3 ) = G (ui,U2,u 3 ) x(ui)fc(u2,u 3 ), , lg , 
H(vi,V2,v 3 ) = G (vi,V2,v 3 ) x(vi)U{v2-,v 3 ). 



Furthermore, by considering a zero-range potential, with 
X = 1, the above coupled equations (17) are reduced 



to the following, where the need of regularization in the 
momentum integrals is explicit: 



K(u%, U3) = 47rr(£ u ) / du' 2 U2 / dx 







f\x 









+ \ j dx> | ( ^ f ni ( w 2) u 2,a;),n 2 (tt2,W3 ) a;')' n3 ( u 2' U3 ' a:; ') J^fn 2 («2>M3,a;'),Il3(M2,W3 I a;') 
+ Go I IIi(w 2 ,«2,a;) ) n4(u2,W3,a;'),n5(w2 ) W3,a;') )%( U4(u' 2 ,u 3 ,x') ,U 5 (u' 2 ,u 3 ,x') 



(19) 



U(v 2 ,v 3 ) = 4ttt(£ v ) / dv' 3 v 3 



J dxGo(^V3,Ue(v2,v' 3 ,x),IlT[v2,v' 3 ,X 



Following an appropriate renormalization scheme for 
zero-range interactions, with the regularization param- 
eters directly associated to observables [2J, in the regu- 
larization of the above coupled equations we introduce 
the short-range three- and four-body scales, as we are 
going to discuss in the next subsection. 



C. Four-boson zero- range bound-state model 

In order to help the understanding of the zero-range 
model in the four-boson formalism, we shortly review the 
main steps to deal with two- and three-boson equations, 
when we have two-body zero-range interactions. The mo- 
mentum space representation of a zero-range interaction 
V(r) — (27r) 3 A <5(r), is characterized by the two-body 
coupling constant A and a separable potential with con- 
stant form factor ((p|x) = !)• The matrix element of 
the two-body transition operator, (p'\t(£)\p) = t(£), is 
obtained by analytical integration of Eq. |T5| ) . The cou- 
pling constant A is fixed by one physical input, e.g., the 
position of the two-body pole at E2: 



A" 1 = / d 3 p 



1 



E 2 -P 2 



cfp 



1 



E-> 



P 



(20) 



We consider bound or virtual two-body energies, when 
E2 = — \E 2 \- Bound in the case the pole is in the up- 
per part of the complex imaginary axis of the momen- 
tum plane (ikb), with corresponding energy in the prin- 



cipal Riemann sheet (E 2 



j); and virtual when 



axis (-ik v ), with the energy in the second Riemann sheet 
(E 2 = As we are going to consider one case or the 

other, both Kb and ny will be labelled as y/B^- There- 
fore, the corresponding scattering lengths are given by 
a = ±l/\fB~2 (+ for bound, and — for virtual), with the 
renormalized two-body <-matrix given by 



t(£) 



1/(2tt 2 ) 



1/(2tt 2 ) 



±l/o • 



(21) 



the pole is in the lower part of the complex imaginary 



The above procedure is enough to render finite the scat- 
tering amplitude, providing the cancelation of the linear 
divergence of the momentum integral. 



1. Subtractive regularization technique 

The subtraction technique used to define the two-body 
scattering amplitude can be generalized to three-boson 
systems. The need for a new parameter beyond a is de- 
manded by the Thomas collapse 48 of the three parti- 
cle ground state, when the two-body interaction range 
ro goes to zero. Such collapse is also associated to 
the Efimov phenomenon by observing that in both the 
cases |a|/f*o - ► 00 [49]. The Thomas collapse is avoided 
at the expense of regularizing the kernel of the three- 
body equation, which is done by introducing a subtrac- 
tion at an energy scale — fi 3 . Together with the two- 
body scattering length, this procedure determines the 
low-energy three-boson properties. The regularization 
of the momentum integration is done by the subtracted 
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form Gq(E) — Go(-fil) in substitution to the free three- 
body Green's function Gq(E). The regularizing energy 
parameter /x| is also called "three-body scale" in view of 
its direct association to a three-body physical observable 
(scale) in the renormalization procedure. 

By adding a fourth particle to the three-body system, 
within the same zero-range two-body interaction, another 
regularization is required in the corresponding formalism, 
due to new terms in the coupled integral equations, not 
directly identified with the three-body kernel. This reg- 
ularization is followed in a similar way as done in the 
three-boson case, by introducing one more regularizing 
parameter ("four-body scale"), fi\, in the FY formalism. 
This new scaling parameter appears in the integrands as- 
sociated to the presence of the fourth particle, in order to 
allow the complete regularization of all the momentum 
integrals. Among the 18 FY components, only the first 
three (i.e. ijk + l,jki + I and kij + I) will fully describe 
the three-body (3B) system (ijk), where the 3B scaling 
parameter /it§ enters in the subtracted form of the free 
Green's function. The 4B scaling parameter /i| enters 
in the subtracted form of the Green's function that are 
present in the remaining 15 components. Therefore, in 
the three components associated to the 3B system, the 
regularization is done by 



G (E) — > G { 3) (E) ee G (E) - G (-^), 



(22) 



other fifteen components: 

G (E) — > G^iE) ee G (E) - G (-/4). (23) 

Summarizing our regularizing approach: The physical 
scales of the four-body problem are the energies of the 
dimer, one reference trimer and one reference tetramer. 
The subtractive procedure adds the regularization pa- 
rameters /U3 and /i4 to the kernel of the FY equations, 
which are correlated to the physical three- and four-body 
scales. The dependence on the subtraction points is elim- 
inated in favor of two observables, which in our case are 
the values of the trimer and tetramer reference energies. 
The elimination of the dependence on the subtraction 
point is possible as evidenced by the existence of corre- 
lations between four-boson s— wave observables. When 
presented in terms of dimensionless quantities, the cor- 
relation between two consecutive tetramer energies is in- 
deed given by a limit cycle. This correlation, presented 
in Ref. |26j and further explored in the present work, is 
in fact consistent with results obtained by other models. 

The present four-body regularization strategy keeps 
fixed the minimal requirement of the three-body prop- 
erties, while the further freedom is used to introduce the 
new s cale . Therefore, in the regularization procedure of 
Eqs. (19), by introducing the required three- and four- 
bod y re gular izat ion parameters as explained above, us- 



with the new scale appearing in the regularization of the 



ing ( 22 1 and ( 23 ) , we obtain the corresponding zero-range 
subtracted FY equations, which are given by: 



K,(u2,u 3 ) = A-kt(£u) I du' 2 u 2 i dx 



G ( Ui(u' 2 ,U2,x),u' 2 ,U3 tC( U 2 ,U 3 



+ \ j dX ' G ^ ) (^^{ u> 2,U2,x),^2(u 2 ,Uz,x'),YIz(u 2 ,U3,x')^k(^ 

+ \ j dX ' G 4) (^l(u 2 ,U 2 ,x),U i (u' 2 ,U S ,x'),n 5 (u 2 ,U3,x')^H^4(u'2,U3,x'),U 5 (^ 

U(v2,V 3 ) = 4tTt(£ v ) J dv' 3 V 3 GJ) 4) (^V3,V2,v' 3 \'h(vi,v' 3 \ 

+ J dxG { 4) ^v 3 ,n 6 (v2,v 3 ,x),n 7 (v 2 ,v 3 ,x)^lC^U(i(v2,v' 3 ,x),Il 7 (v2,v' 3 ,x)^j 



(24) 



The physical picture behind the regularization leading 
to (24) can be described as follows. The three-body scale 



parameterizes the short-range physics in the virtual prop- 
agation of the interacting three-boson subsystem within 
the four-body system. The four-body scale parameter- 
izes the short-range physics beyond the three-body one, 
explored in the four-boson virtual-state propagation be- 
tween two different fully interacting three-body or dis- 
joint two-body clusters. It is worthwhile to stress that 



the first term in the right hand side of the equation for 
1C(u2,us) is the three-boson bound-state equation, when 
all the other terms are dropped out. So, it is expected 
that such term should carry the 3B scale, while one could 
think that the other terms do not require regularization. 
However, as shown in [T5], still all the other terms need to 
be regularized, in order to avoid the 4B ground state col- 
lapse. Although the same regularizing parameter could 
also be used in all the terms, the freedom to choose it 
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differently is at our hands, and the constraint of intro- 
ducing the trimer properties does not require any further 
assumption, but the given form of the first term in the 
right-hand side of the equation for £(^2,1*3). The phys- 
ical conditions in cold atom traps that allow us to ex- 
plore a short-range four-boson independent scale, when 
approaching the Feshbach resonance, are discussed in 
Ref. HI]. 



III. BINDING ENERGIES OF UNIVERSAL 
TRIMERS AND TETRAMERS 

In this section, before considering our numerical re- 
sults obtained for the FY bound-state equations ( 24 ) , we 



provide a brief discussion on scaling functions, starting 
with the three-boson case (next subsection) and mov- 
ing to tetramer systems, in subsection III B| After we 
set up the concept of a scaling function, we will present 
detailed numerical results for the tetramer binding ener- 
gies, with a comparison with other recent available calcu- 
lations. The Lanczos-type procedure for solution of the 
coupled FY equations is shown in Appendix [C] where we 
also give some details on stability and convergence of our 
numerical approach. 



A. Trimers in the scaling limit 

The three-body system is sensitive to the physics at 
short ranges, which is parameterized by /13. After this 
scale was recognized [2] , it was shown that the change of 
the three-body scale in respect to the two-body one (e.g., 
1/et) can be clearly revealed by expressing the energies 
of the Efimov trimers in a single curve, which defines 
a scaling function T^ N \ It allows to build the energies 
of all the sequence of weakly bound trimers. Thus, the 
energies of successive trimers are correlated by [2] 



\ 



B 



{N+l) 



B. 



(AO 




(25) 



written in term s of dimensionless quantities. For conve- 
nience, in (25), we define B 2 



systems (p 



B 2 for bound two-body 



us sign, or a > 0), and B2 = for virtual 
states (minus sign, or a < 0). 

Few cycles are enough to reach a universal function 
independent on N; i.e., in the limit N — > 00 it reaches 
a renormalization-group invariant limit cycle pQ. At the 



critical values, B 



(N+l) 



B 



(N+l) _ 



= B2 (for bound two-body) and 



(for virtual two-body), the scaling function 



(25) vanishes: F 3 
the 



(AT) 



± 



= 



In this limit, 



solution for a bound two-body system (+ sign) is 
5 B 2 ; and the solutio 
two-body system (— sign) is given by B^ 



given by £3 — 6.925 B 2 \ and the solution for a virtual 



} 3 



where B 2 is the virtual state energy [2J [50] . We observe 



that the scaling law (25) is one among many possible 
model independent correlations between three-body ob- 
servables for short-ranged interactions. We should also 
point out that range effects can become relevant as |a|/ro 
decreases. This is evidenced by the finite-range numer- 
ical results from other authors, presented in the scaling 
plot displayed in Ref. [2J. (On range effects and univer- 
sal properties of three-body systems, see also Ref. [51] ■) 
Therefore, range effects will affect Eq. (25 1 and the above 



(AT) 



critical values of B^ 

Evidence of Efimov cycles, with a period of 22.7, in 
the values of the two-body scattering length, a, at the 
peak of the three-body recombination, were observed in 
a cold-atom experiment by the Innsbruck group (25], as 
well as in other recent experiments reported in [9j . These 
results gave a response to a quite old discussion on the 
possible realization of the Efimov effect in real physical 
systems. The Efimov effect was formulated in the nu- 
clear physics context, when considering the solutions of 
the three-body Faddeev equations in the limit of infinite 
large two-body scattering length. Since then, the possible 
existence of Efimov states in nature have been considered 
in several works, where the most promising one was the 
calculations done by Cornelius and Glockle for a system 
with three atoms of 4 He [52] . The revival of this matter 
came recently with the advances in cold-atom laborato- 
ries, where by using Feshbach resonance techniques it is 
possible to vary the two-body scattering length from zero 
to infinite (positive or negative) values. 



B. Tetramers in the scaling limit 

Following a complete analogy with the three-body 
case, we now show evidence for the existence of a limit 
cycle for the binding energies of successive tetramers. 
Therefore, by considering a shallow dimcr (bound or vir- 
tual), we introduce a scaling function for the energies of 
two successive tetramers between two trimers, as 



\ 



B 



(N+l) 



B, 



B 



(AO 




(26) 



For the present purpose, the calculations are restricted 
to bound tetramers below the ground trimer, although 
the results are valid for tetramer attached to different 
trimers [26 . 



The scaling function ( 26 ) has the dependence on a four- 



body scale independent of the three- and two-body ones. 
Due to the flexibility of our model it is transparent how to 
obtain it. The FY equations (24) are solved for /i 3 fixed 
while fi4 is varied. In the unitary limit, the solutions 



of (24) depends only on the ratio /U4//X3, with observ- 
ables given in units of /13. Eliminating the dependence 
of the tetramer energies on the ratio ^4/^3 by writing 
the energy of the (N + 1)— th state as a function of the 



energy of the N— th tetramer, the scaling function (26) 
is constructed, where the dependence on /13 is removed 
in favor of B 3 . As no other scales are present in (24 1, a 
scaling function like (26) seems possible. It remains to 



find whether Eq. (26) is independent on N, which will be 
discussed in the next subsection. 

We observe that, tetramers below a generic trimer 
are constrained to be between the two successive Efi- 
mov trimers. However, below the ground state trimer 
the tetramer collapses as 00 |15) . as it is exem- 

plified in Table [I] (A) , for scale ratios above 200 in the 
unitary limit. 



C. Numerical results close to the unitary limit 

The numerical calculations with the renormalized zero- 
range model requires the definition of the relevant scales. 
They correspond to the two-, three- and four-body scales 
given by a -1 , /i3 and 114, and obviously for fixed scatter- 
ing length and trimer scale one can move ^,4 to investigate 
its effect on the spectrum of tetramer ground and excited 
states. 

In Table ||J we have listed our numerical results for 
tetramer ground and excited state binding energies. In 
part (A) the results at unitary limit B2 — are given 
for different scale ratios from ^4/^3 = 1 to 400. Accord- 
ing to the obtained results for tetramer binding energies, 
/X4//Z3 ~ 1.6 is the threshold for the first tetramer excited 
state, 1^4/^3 — 21 is the threshold for second tetramer ex- 
cited state and close to /U4//U3 ~ 240 the third tetramer 
excited state should appear. At these critical values of 



scale ratios where B 



(JV+l) 



function ( 26 ) vanishes 



= B 3 with B 2 = 0, the scaling 




0=0 



with the solutions approaching B^P C ~ 4.6 B%, for each 
cycle, i.e., N =0 and 1, as we have calculated up to three 
tetramers below the trimer. The limit cycle is approached 
quite fast, as verified in the case of three-boson [5] . In the 
critical condition (27) we have one tetramer at the atom- 



(27) 



trimer scattering threshold, which allows for a resonant 
atom-trimcr relaxation in trapped cold-atom gases at the 
Feshbach resonance. 

In Table [j](B) and (C), we show the numerical results 
for posi tive an d negative scattering lengths, respectively, 
with \[B~?JBj, = ±V0.02. From these results it is clear 
that, by comparing with the case of infinite scattering 
length and for the same scale ratios, the tetramer bind- 
ing energies for positive and negative scattering lengths 
have smaller and larger values, respectively. The critical 
condition for atom-trimer resonant relaxation near the 
Feshbach resonance, has to be corrected as 



(28) 




TABLE I: Binding energies of ground and excited tetramer 
states, for different four-body scales and for (A) B2 = 0; (B) 
bound dimer, with B2 = 0.02 B3; and (C) virtual dimer, with 
\fB~2~ = — V0.02 B3. In (A), we verify that a third excited 
tetramer emerges for 1x4/113 w 240. In (B) and (C), for non- 
zero two-body binding, we have only presented results for the 
ground and first excited state binding energies. 



(A) 



II A / lit 


B^/B-i 

J -'4 / 


J - / 4 / •-> 


- 1 B, (2) IB% - 1 


1 


3.10 






1.6 


4.70 


7.10x10 


-4 


5 


12.5 


0.531 




10 


24.6 


1.44 




21 


63.5 


3.62 


3.20xi0" 4 


40 


184 


7.65 


0.203 


50 


275 


9.80 


0.365 


70 


520 


12.9 


0.629 


100 


1.04xl0 3 


20.5 


1.17 


200 


4.06 xlO 3 


50.8 


2.86 


300 


9.11xl0 3 


102 


4.53 


400 


1.62xl0 4 


153 


6.28 








1 


2.66 






1.76 


4.24 


9.8xl0 -4 




5 


10.0 


0.421 


(B) 


20 


45.9 


2.77 


40 
80 


139 
506 


6.10 
13.0 




200 2.86xl0 3 


39.5 




300 6 


OOxlO 3 


69.3 




400 9.81xl0 3 


104 



(C) 





</B 3 




1 


3.62 




1.7 


5.91 


0.014 


5 


15.4 


0.658 


20 


74.8 


4.18 


40 


236 


9.46 


80 


873 


20.6 


200 


5.02xl0 3 


64.5 


300 


1.06 xlO 4 


115 


400 


1.73xl0 4 


174 



a = ±l/-*/B~2 does not vanishes. By performing a lin- 
ear expansion around a -1 = 0, from the results given in 
Table |(B) and (C), one has 



B 



(JV) 



4.6 B 3 



1 - 0.8 {ayflh)- 1 



(29) 



The dependence of the tetramer energies with the ratio 
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FIG. 2: (Color on-line) The first three tetramer energy levels 
at unitary limit, in units of the trimer ground state, are shown 
as functions of ((14/ (13) 2 , where /Lt4 and /13 are, respectively, 
the four- and three-body regularizing scaling parameters. 



/X4//X3 in the unitary limit is presented in Fig. [2] by using 
results of Table |l] (A). The collapse of the tetramer states 
by increasing the short-range four-body momentum scale 
in respect to the three-body one is seen. As for 1/a = 
0, only two scales define the tetramer state, the four- 
boson binding increases in respect to the trimer one, by 
either increasing /i 4 or decreasing /i 3 . The three tetramer 
states for N = 0, 1 and 2 show an increase of the binding 
energies roughly as A similar scaling appears for the 
trimer binding energy with The ratio of B4/B3 gives 
the slope of the Tjon line and it is not a constant as 
thought before. As a matter of fact it can be considered 
as a parameter measuring the four-body scale. 

In order to compare directly the numerical results for 
the three- and four-body limit cycle, we present the Fig. [3] 
with two panels, (a) and (b). The three-body Efimov 
states are shown in (b), with the energies given in units 
of the trimer scaling parameter The corresponding 
four-body Efimov-like states are given in (a), when con- 
sidering a ground-state trimer in the unitary limit. The 

dependence of ZJ4 with B3 are given with the energies 
in units of the tetramer energy parameter for N =0, 
1 and 2 (with fixed trimer parameter /i 3 = 1). We il- 
lustrate the main qualitative phenomenon of tetramers 
emerging from the atom plus trimer threshold: by de- 
creasing £>3, an increasing number of tetramers become 
bound. A similar phenomenon is shown in the right panel 
for trimers in terms of a dimer energy, which is the Efimov 
effect. In our calculations, we first consider tetramers be- 
low a trimer in the ground state, such that the tetramer 
spectrum is not limited from below and can collapse as 
we increase the four-body scale in relation to the three- 




10"" 1(P 10^ 1 

B 2 / ii 2 3 

FIG. 3: (Color on-line) In (a) we have the tetramer energy 
levels N — 0,1,2, as functions of the trimer ground-state 
energy for B2 = 0. For comparison, in (b), we have the first 
Efimov trimer levels as functions of B2 ■ The energy units are 
pf, in (a) (where /i§ = 1); and /^|, in (b). In both cases, the 
threshold is given by the diagonal line. 



body one. In this case, the number of excited tetramer 
levels can increase with no limit. However, for a given 
general excited trimer, the applicability of our results are 
restricted to a region where the trimer energy varies by 
a factor of about 515 (the Efimov ratio in the unitary 
limit). Therefore, between two successive trimer states 
at most three tetramers can be found. 

The results for the four-boson binding energies, plot- 
ted in Fig. [3^a), exhibit a limit cycle, which expresses 
the universal behavior of the energies with a moving 
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T , , 1 1 , 1 1 , 1 1 , 1 1 1 1 1 , 1 1 , 1 r 




FIG. 4: (Color on-line) The four- and three-body scaling func- 
tions are represented, respectively, by M = 4 and M = 3. The 
energies are rescaled in terms of Bj, where -/V represents the 
corresponding energy level in both the cases (N = is the 
ground state). The behavior of excited trimer (M = 3) and 
tetramer (M = 4) energies are shown as the corresponding 
subsystem energy ( B2 or B3 ) is varying. In the case of 
M — 4, B2 is fixed to zero and the dashed part is only reach- 
able by tetramers below the ground-state trimer. 

four-body scale parameter. The curves shown in this 
figure reduces to a single one, when they are plotted as 
the correlation between successive tetramer energies, as 
presented in Fig. [4] for M = 4, where we consider the 
scaling plot in the unitary limit. The scaling function 

(J B^/B^; O^j was built considering up to the 

third excited state. As in the case of trimer Efimov cy- 



cles, the numerical results for tetramers also present a 
very rapid convergence towards a four-boson limit cycle. 

For comparison, we also present ^Jb~2~/B^j , 

which has the value of 1/22.7 when B2 vanishes. The 
three-body scaling function, identified inside the plot by 
M = 3, was derived in Ref. [2J, being compared directly 
with the corresponding four-body scaling function, iden- 
tified by M = 4. The dashed part of the tetramer scaling 

curve presented in Fig. j^J where \J B^/B^ 1 ^ < 1/22.7, is 
not accessible for the excited tetramer energies belonging 
to the tetramer spectrum that is obtained for a given ex- 
cited trimer. This restriction implies that a meaningful 
infinite number of tetramer levels is only possible when 
considering the trimer in the ground-state level. 



Let us consider the unitary limit (when we have the 
geometrically spaced trimer levels), to resume the global 
picture suggested by our model. By increasing the 
tetramer scale in respect to the trimer one, as a new 
excited tetramer emerges from the threshold, the less ex- 
cited ones will slide down, moving in direction of the next 
deeply bound trimer. Actually these tetramers, between 
two trimers, are resonances decaying to an atom and a 
trimer. They will dive into the complex plane through 
the atom-trimer scattering cut with their width increas- 
ing as the four-body scale is incremented. This qualita- 
tive discussion deserves further investigations from the 
point of view of the complex analytical structure of the 
four-body scattering equations. Our picture suggests 
that, in the unitary limit, at most three tetramer res- 
onances lie between two successive trimers. 



In Fig. [5] to verify the universality of the four-boson scal- 
ing function ( 26 1 and for comparison, we include results 
obtained in other recent calculations. In Table ILT1 we have 
listed the available results of other authors, for binding 
energy ratios of ground and excited tetramers close to the 
unitary limit. Within those results we should also point 
out that the ones given in Refs. [TTlQjl] are for more than 
one trimer levels. In our results with zero-range interac- 
tion, we include not only the exact limit B2 = 0, but also 
non-zero two-body energies (bound or virtual) , as shown 
in Tables |l| (A), (B) and (C). The sensibility of the scal- 
ing function with variation of the two-body conditions is 
shown when non-zero two-body energies are used. As it 
is shown, the exact unitary limit scaling plot is shifted 
to the right-hand side (left-hand side) when considering 
non-zero two-body bound (virtual) state energies. 

The results of available calculations for tetramer ener- 
gies plotted in the way we are suggesting put in evidence 
the effect of the four-body scale. Consistent with our 
findings, they slide along the universal correlation shown 



in Fig. [5] The sensitivity of our defined scaling plot, in 
respect to variations of a two-body observable, as the 
scattering length or range effects, is also consistent with 
the results obtained by other authors. 

In order to verify the sensitivity of the results (ob- 
tained for B2 = 0) with small changes in the dimer ener- 
gies, we calculate the corresponding scaling functions for 
a few cases with B2 7^ 0. In the given results, we have 
three cases for bound two-body systems, where B2/B3 = 
0.02, 0.0044, 0.002, and on e cas e for a virtual two-body 
state, with ^B 2 /B z = -\ZO02. In Fig. |HJ we have also 
included the Hammer and Platter results for both vir- 
tual and real dimer close to unitary limit (See Fig. 2 
of Ref. [H]). Despite of the observed difference in the 
position of both results in our plot, it is remarkable that 
they both fit nicely in our scaling curve, meaning that the 
four-boson model of Ref. [TH] is sensitive to a four-body 
scaling parameter, carried out implicitly by the momen- 
tum cut-offs used in their calculations. 

Although we calculate tetramer energies below the 
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FIG. 5: (Color on-line) Our results for the four-boson scaling function, reported in [26], are displayed in the limited region 
where more recent results are available. The plots are shown for B2 = (solid-black curve) and in four cases with B2 7^ [three 
cases with bound dimers, B2/B3 — 0.02, 0.0044, 0.002, and one case with virtual dimer, ^/B^jBz = — \/0.02], as indicated 
inside the frame. In the case of B2 = 0, this figure refers to the right-hand-side corner of Fig. [4] for M = 4. The results were 
explicitly verified numerically for N = and N = 1, considering the energies B^ ' 1 ' 2 , in case B2 = 0; and for N = in the 
other cases where B2 7^ 0. The symbols refer to other model calculations (as indicated) near the unitary limit. 



trimer ground state, such that the dashed part of Fig. [4] 
is also verified, the results for tetramers attached to any 
other Efimov state, does not change the present conclu- 
sions on the existence of a proper four-body scale or its 



universal manifestation through the scaling function ( 26 1 , 



calculated with the zero-range model regularized within 
our scheme. This fact is evidenced in Fig. [5] Distinct 
short-ranged interaction models [TBI H7l IT9H2T| show that 
the energies of successive tetramers linked to the tower 
of Efimov states scale according to the plot, and thus 
verifying the universality of the four-boson limit cycle. 



IV. STRUCTURE OF UNIVERSAL 
TETRAMERS 

The effect of the four-body scale on the structure of 
tetramers at the unitary limit is detailed in this section. 
We show results for the FY components in momentum 
space for different scale ratios as well as for momentum 



probability densities. A close inspection to the reduced 
FY components, K, and %, will show how 114 manifests 
through their dependence on the different Jacobi mo- 
menta. An analogous study will be presented for the 
momentum probability densities. 



A. Reduced Faddeev-Yakubovsky Components 

The reduced FY components IC(u2,Us) and / H(v2,V3), 
are shown to spread out up to momentum of the order 
1, due to the relevance of the four-body momentum scale 
/J4 to regularize the kernel of the set of FY equations 
(24) at short distances. Both components have narrow 



peaks appearing at small momenta, which are even more 
pronounced for higher scale ratios and more ex- 

cited states. They follow a typical tetramer momentum 

scale ( 
and [7] 



B[ N ^ /Hi). These features are evident in Figs. 
For a given ratio ^4/^3, the ground state wave 
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TABLE II: Binding energy ratio of tetramer ground and 
excited energies to a trimer level, ^B^/ife, B^/B 3 — 1 
where the results of Refs. I17l 119] refer to different trimers. 



Ref.[53] 


Ref.[T4] 


Ref.HI] 


Ref.[E] 


Ref.[20] 


[4.46,0.06] 


[4.075,0.003] 


[4.41,0.01] 


[5.0,0.01] 


[4.55,0.003] 



Ref.rm 



Ref.Ojj] 



[4.58, 0.01] 
[5.88, 0.10] 
[4.58, 0.006] 
[5.58, 0.03] 
[4.55, 0.001] 



[4.5175, 0.00106] 
[4.6041, 0.00217] 
[4.6104, 0.00227] 
[4.6108, 0.00228] 
[4.9929, 0.00997] 
[4.6114,0.00228] 



function is more expanded in momentum space in com- 
parison to excited states. The K— and i7-channels, and 
consequently the total wave function for smaller tetramer 
bindings, should extend to larger distances as compared 
to the corresponding ones with larger binding. 

The peaks at small momentum are due to the nearest 
trimer pole in the three-boson interacting resolvent, ap- 
pearing in the A"— channel equation (24 1. This leads to 



a dominance of K, over H for small momentum when the 
tetramer energy approaches B 3 . 

In Figs. [6] and [j] we detail the sections of JC(u 2 ,u 3 ) 
and H(v 2 , Va) with one of the Jacobi momenta being zero 
for scale ratios of 5, 50 and 200. The interpretation of 
the Jacobi momentum in respect to the relevant scales 
is necessary to proceed in the analysis of these plots. 
The variable u 2 is the relative momentum of particle 3 in 
respect to the center of mass of the pair (1,2), which be- 
longs to a trimer configuration (see Fig. [I]). Note that in 
our derivation of the FY equation for the zero range po- 
tential, we have chosen to factorize out the dependence 
of the relative momentum of the pair (1,2) in the defi- 



nition of the reduced amplitudes, as given by Eq. (18). 
Naively, by invoking kinematical reasons, it is reasonable 
to expect that, within the tetramer, the average distance 
between two bosons is larger than the average distance 
between the center-of-mass (cm.) of the virtual pair to 
the third particle, which is larger than the average dis- 
tance of the virtual trimer to the fourth particle. The 
ordering effect comes because the clusters being heavier 
should have their cm. closer to the cm. of the tetramer. 
This size ordering has a correspondence in respect to the 
momentum variables. In K. (112,1*3), the momentum of 
the fourth particle in respect to the virtual trimer, U3, 
should explore larger momentum regions than u 2 , which 
is the relative momentum of the third particle in respect 
to a pair. An analogous reasoning suggests that the mo- 
mentum dependence of T-L(v 2 , V3) is such that v 2 explores 
larger momentum regions than V3 (the relative momen- 
tum of two bosons). These qualitative properties are ver- 
ified in Figs. [6] and [7] 
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FIG. 6: (Color on-line) The Yakubovsky components JC 
and T-i, as functions of the Jacobi momenta for scale ratio 
^4/^3 = 5, when only one four-body excited state exists. 
In the frames (a) and (b) the components are shown as a 
function of 112 , V2 where u 3 , 113 = and in the frames (c) 
and (d) are shown as function of 113,1)3 where U2,V2 = 0. 
In (a) and (c) we have results for the ground-state level; 
and in (b) and (d) for the first excited-state level. The 
normalization is such that J °° du2 u| J °° du 3 u§ K. 2 (u2, 113) + 
f °° dv 2 v\ f °° dv 3 «| K 2 (V2, v 3 ) = 1. 



The plots in Fig. [6] show the reduced FY components 
for the ground, frames (a) and (c), and excited state, 
frame (b) and (d), for ^4/^3 = 5. The high momen- 
tum tails of H(v 2 ,0) and /C(0,U3) in respect to 1-1(0, V3) 
and K,(v,2, 0) are visible by comparing the frames (a) and 
(c) in the case of the ground state, and (b) and (d) for 
the excited state. By inspecting the pair of frames (a) 
against (b) ; and (c) against (d) , we also note that all the 
FY components for the excited state are more concen- 
trated at lower momentum than the corresponding ones 
for the ground state. Although the binding energy ratio 
/B3 for the excited state is close to 1, the high mo- 
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FIG. 7: (Color on-line) Following the caption of Fig. [6] here 
we have the same quantities for n^/fis = 50. In the frames 
(a) and (d) we have the ground state; in (b) and (e), the first 
excited state; and in (c) and (f), the second excited state. 

mentum tail bringing the four-body scale is present, giv- 
ing to this state the possibility to move as (14 is changed 
in respect to /x 3 . These qualitative features are supported 
by the results we have obtained for larger scale ratio, such 
as ha/ ^3 =50 and 200, where three tetramer states are 
possible. This behavior is clearly shown in Fig. [7] for 

In all the cases we analyzed, the dependence of U3 in 
K{u2,U3) and V2 in 7l{v2,V3) show that these two vari- 
ables explore larger momentum regions than ui and V3, 
respectively. As we have shown /C(u 2 ,0) and 7^(0, i> 3 ) 
have a less prominent tail at large momentum than 



£(0,1(3) and H(v2,0)- If, for any reason, the high mo- 
mentum dependence that appears in the K— and H— 
channel are dropped out, the dependence on the four- 
body scale will be minimized (may be even completely 
removed!). Therefore, unreasonable selection of cut-off 
values in the mapping of momentum variables, which can 
happen in particular when a fixed cut-off is being used 
for increasing values of the scale ratios, can lead to con- 
vergence in the four-body binding energies rather than 
collapse. For the discretization of momentum variables 
one should not only consider large enough cut-off val- 
ues, consistent with the used four-body scale, but also 
consider a reasonable number of mesh points in the in- 
terval near zero-momentum. Since the iteration of the 



coupled Eqs. (24) requires a very large number of multi- 
dimensional interpolations on the Yakubovsky compo- 
nents, we have used Cubic-Hermite Splines to reach high 
computational accuracy. 

From the above analysis, represented in Figs. [6] and [7j 
we conclude that, against the general belief stating that 
in four-body atomic and nuclear calculations, with model 
potentials in the FY scheme, the K— channel is always 
dominant, we show some cases where the H— channel is 
dominant even at low momentum. This is shown in par- 
ticular for high scale ratios /Z4//U3. 



B. Momentum probability densities 

The tetramer wave function ^(iti, 1*2,^3) fully sym- 
metrized is built from Eq. ^ using the definitions ( 18 ) 
of the K and H components in terms of /C(it2,M3) and 
"H(v2, V3). In order to simplify our analysis of the wave 
function, and to obtain an insight on how the momentum 
is shared among the Jacobi coordinates, we introduce the 
momentum probability densities n(ui) as: 



n(ui) = iq J dujU 2 j j duku\ ^ {u\, 112, 1*3) ;(30) 

f° 

Jo 



JO 

dui n(ui) = 1, 



with (i,j,k) = (1,2,3), (2,3,1) and (3,1,2). For reference 
we use the variables Ui of the K configuration, where 
Mi is the relative momentum of a pair, u 2 the relative 
momentum of the third boson to the CM. of the pair, 
and U3 the relative momentum of the fourth boson to 
the CM. of the other three (see Fig. [T]). As we have 
discussed in sect. |IV A| the average values of the mo- 
mentum variables arc ordered as (ui) < (1*2) < (1*3), 
which is reflected in momentum densities, independent 
of the degree of excitation of the tetramer. 

The momentum probability densities n(ui), n(u2) and 
n(u3) for scale ratio ^4/^3 = 50, for ground, first and 
second excited tetramer states are presented in Fig. [8] 
In the left panel we show our results for the probabil- 
ity densities in units of /Z4. Clearly for all three mo- 
mentum densities, the second excited tetramer leads to 
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a very sharp peak close to zero momentum, whereas for 
the first excited and ground tetramers the peak is more 
wide, which is evident as the size is inversely correlated 
to the tetramer binding energy. For the ground tetramer 
state the FY components and consequently the 4B to- 
tal wave function are expanded to higher momentum re- 
gion and for excited tetramers they are more restricted 
to zero momentum region. In the right panel of Fig. [8] 
we present the probability densities in momentum units 

= "J B±*\ for ^4/^/3 = 50. By considering this scal- 
ing factor corresponding to the energy of each tetramer 
state, which from Table [T] varies two order of magnitude, 
the densities are amazingly close. The tetramer energy 
can be considered as the physical scale, which correlates 
other observables. 

One feature of the universal properties of the four-body 
wave function is presented in Fig. [9] , where we have com- 
pared the momentum distribution functions n(ui) of the 
two shallowest tetramers, for two very different values 
of ^4/^3 equal to 5 and 50. The functions n{ui) were 
rescaled to the same maximum value to make transpar- 
ent their universal form. The comparison between them 
clearly confirms a universal dependence of the momen- 
tum distributions on the tetramer binding energy. Note 
that, these quantities are plotted in terms of dimension- 
less quantities in natural units of the length associated 
with the tetramer energy. 

In Fig. [9] we have shown the renormalized form of the 
rescaled distribution functions, considering two different 
scale ratios ^ = 5 and 50, whereas the distributions are 

fJ-3 ' 

renormalized to the same maximum value. As we have 
shown the distribution functions have universal shape, 
independent of the value of scale ratio. We should also 
add a remark that, in order to simplify the presentation of 
Fig. [9j we have included only the momentum distribution 
functions for the two more excited states in case of — = 

fJ-3 

50, whereas the same universal shape is also obtained for 
the momentum distribution functions of the ground state 

(see Fig. pj. The ground state energy ratio is B^ j B 3 = 

275 (see Table|l|, which means that is about 2. 7/13, 



while for the two excited states, B^ 1 and B^' is much 
smaller than thus the three-body regularization scale 
is far more important to the ground state than for the 
excited states, which approach an universal form for the 
momentum distribution functions. 



? (2) 



ators O with the corresponding FY amplitudes denoted 
by Kq and Hq (see Appendix [B]) . Computing Kq and 
Hq for different O, the 4B total wave function can be 



obtained by using Eq. ( Bl I or (B2| 



For our purposes of presenting graphically the total 
wave function we will show plots for the angle aver- 
aged wave function, where the dependence on the rel- 
ative angles between Jacobi momentum vii (i = (1, 2, 3)) 
is integrated. Our notation is such that we still use 
$?(ui, 1*2, M3) for the angle averaged wave function. In- 
stead of the usual normalization, we consider 



<*!*> 



du\u\ 



du 3 ul * 2 (ui,m 2 ,u 3 ) 
(31) 



In Figs. [T0p2| we present our numerical results for 
the angle averaged total wave function with scale ratio 
= 50, for the ground, first and second excited 
states. The wave functions are shown as functions of two 
Jacobi momenta, when the third one is chosen to be zero. 
As we did before, the magnitudes of Jacobi momenta 

are re-scaled by the momentum factor (5^ = \J~B^ 
(N = 0,1,2). In view of the normalization condition 
(31 1, this momentum scaling implies in a rescaling of the 



magnitude of the corresponding total wave function by 
the factor of /3^ 2 . 

The universal form of the 4B wave function is evi- 
denced in Figs. T0][T2"] by a proper rescaling with the 
four-body energy. Correspondingly, the contour plots are 
also presented as functions of two non-vanishing Jacobi 
momenta, given in units of /3/v- Also the momentum den- 
sity distributions, for two excited states, present a univer- 
sal form, i.e., independent on the binding energies, when 
the momentum is measured in units of j3jq, as shown in 
Fig. [8] Note that each of the choices of the vanishing mo- 
mentum implies in a very long wavelength limit for the 
fourth particle, indicating the dominance of a three-body 
cluster configuration within the four-body system. This 
configuration carries the four-body scale information, as 
clearly shown by the plots for the wave-function depen- 
dence on the two non-vanishing momentum in units of 



CONCLUSIONS AND PERSPECTIVES 



C. Tetramer wave function 

The calculation of the FY amplitudes K and H, given 
by Eq. ( 18 1 for form factor unity, with the reduced am- 



plitudes JC and H from the numerical solution of the 
zero-range FY integral equations (24), allows to build 



the four-boson wave function fully symmetric by boson 
exchange. The necessary permutations to reconstruct the 
4B wave function from the computed K(u\, U2, M3) and 
H{v\, 1*2, ^3) are given by introducing permutation oper- 



The momentum-space structure of the FY components 
of weakly-bound tetramers was thoroughly analyzed at 
the unitary limit using a renormalized zero-range two- 
body interaction, with three and four-body scales. Our 
regularization scheme allows to introduce independently 
these two scales providing general framework to inves- 
tigate the universal properties of tetramers, which has 
extended tails in the classically forbidden region. The 
step to renormalized results for the observables is found 
by introducing scaling functions, which are written only 
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in terms of dimensionless physical quantities. These uni- 
versal scaling functions achieve numerically a limit cycle. 

The universal scaling function correlating two succes- 
sive tetramer energies attached to one trimer comes from 
the sensitivity of the four-boson system to a short-range 
four-body scale. Each excited N— th tetramer energy 
moves as the short-range four-body scale changes 
while the trimer properties are kept fixed. We suggest 
that the ratio of B4/B3, which gives the slope of the 
Tjon line, and it is not a constant as thought before, as a 
parameter measuring the four-body scale. At the unitary 
limit, the successive (N + 1)— th tetramer state emerges 
from the atom plus trimer threshold for a universal ratio 
b[ N) /B 3 = b{ N) /b[ N+1) ~ 4.6, which does not depend 
on N . The atom-trimer relaxation resonates when the 
tetramer hits the scattering threshold. 

We also find out that other model results obtained at 
the unitary limit or close to it [El [I7J IT9H21"] , are quite 
consistent with our four-boson scaling plot, giving con- 
fidence on the universality of our proposed scaling func- 
tion. As evidenced by our results, the independent be- 
havior of the four-body scale is verified in particular when 
a universal excited four-boson state pumps out from the 
atom-trimer threshold as the four-boson parameter is 
driven to short distances or to the ultraviolet momen- 
tum region. In the case that both scales are similar, we 
confirm that our scaling approach is consistent with the 
results of other model calculations. We note that the re- 
sults obtained by other groups appear near the threshold 
region due to the model assumptions that have been con- 
sidered, which are much too restrictive when compared 
to the flexibility of our model and need to be relaxed to 
allow a wider variation between the four- and three-body 
properties. The suggested scaling plot offers a model 
independent way to view the relation between successive 
tetramer states, which can be verified by experiments ex- 
ploring two-body scattering lengths very close to the Fes- 
hbach resonance, where multi-boson forces arc expected 
to be active [57]. A four-boson short-range interaction 
can drive the four-body scale independently from the 
three-body one. Short-range three-body forces accom- 
panied by different off-shell behaviors of the two-body 
interaction, which kept the low-energy two- and three- 
body properties unaltered, could provide another possi- 
bility to move the four-body scale in practice. 

The attent reader should realize the power of the four- 
body scaling plot in pinning down the independent be- 
havior of a four-body scale (for a fixed three-body one), 
in a similar way as the three-body scaling plot, given in 
Ref. [5] , resumes all the Efimov states obtained by chang- 
ing the three-body scale. Next, we should emphasize 
that the scaling curve relating two successive tetramer 
energies is valid not only for the stable tetramers at- 
tached to the ground-state trimer, but also for the un- 
stable tetramers between two successive trimers. For 
such conclusion we have also performed resonance cal- 
culations of tetramers attached to excited trimers, which 
are under detailed numerical analysis to be included in 



a future work. We avoid to present such results in the 
present work mainly because we have considered other in- 
dependent and precise calculations supporting our claim 
of model independence. In this regard, the recent calcula- 
tions performed by Deltuva [TH] are consistent with our 
results and give confidence on the model independence 
of the scaling plot relating the energies of two successive 
tetramers between two successive Efimov trimers in the 
unitary limit. As we have shown, the effect of scatter- 
ing length variations is noticeable and shows a universal 
pattern in agreement with calculations shown in |16j , per- 
formed in the context of effective field theory. However, 
range corrections, which were not yet explored within a 
general study considering the trimer and tetramer inter- 
woven spectra, will certainly have their importance in ex- 
plaining how results near the unitary limit deviate from 
the calculated scaling function for a zero-range force. In 
this respect, range corrections, linear in the two-body ef- 
fective range, have been already addressed within an ef- 
fective theory treatment of the three-boson problem with 
subtracted momentum-space integral equations |57j . 

The scaling plot representation exhibiting the depen- 
dence on the four-body scale, which is also being con- 
firmed by other models, shows that the model inde- 
pendence is not fortuitous. The reason for that comes 
from the fact that the wave function of tetramers are 
largely dominated by configurations where the bosons 
are outside the potential range in the classically forbid- 
den region and depend just on two scales at the unitary 
limit. The sensitivity to the short-range scales appears in 
both the K— and H— channel of the FY decomposition, 
which present high momentum tails for any degree of the 
tetramer excitation. We also found that the H— channel 
is favored over K— channel at low momentum when the 
four-body momentum scale largely overcomes the three- 
body one. 

The universal form of the wave function is putted for- 
ward at the unitary limit by a scaling plot where the 
probability momentum densities for different Jacobi mo- 
menta are shown in units such that the tetramer binding 
energy is one. We found that the shapes are independent 
on the scale ratio and excitation. It only depends on the 
chosen Jacobi momenta. A simple scaling rule giving the 
ordering (ui) < (^2) < (^3), where u\ is the relative mo- 
mentum of a pair, w 2 is the relative momentum of the 
third particle to the pair, and U3 is the relative momen- 
tum of the fourth particle in respect to the three-body 
subsystem. 

Our calculations of tetramer binding energies, with 
corresponding structure of momentum probability den- 
sity and wave functions, are supplying strong numeri- 
cal evidences that universal tetramers can arise from the 
trimer threshold as the four-body scale is moved. These 
results suggest the possibility of resonant atom-trimer 
recombination process near the unitary limit, i.e., when 
a at ±00. It is worth noting that the problem of a new 
scale in the tetramer properties has a long history, which 
was never clarified before in the physics community to 
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the extend we are doing now. The four-body scale gives 
rise to new scalings of tetramer observables which are not 
determined only by the two and three-body properties. 
We present detailed and accurate numerical results for 
the binding and the structure of tetramer states to sup- 
port our claim, given conveniently within the framework 
of scaling functions. As a matter of fact, even results ob- 
tained near the unitary limit, by authors that are claim- 
ing that no four-body scale is necessary, are consistent 
with the universal scaling plot for the tetramer binding 
energies. 

In our numerical results of four-boson bound states, we 
first consider the exact unitary limit (1/a = 0), and we 
next consider deviations of this limit to bound (a > 0) 
and virtual (a < 0) two-body branches. As shown in 
Figs. 3][5 universal tetramers can hit the atom-trimer 
threshold leading to a resonant relaxation in this chan- 
nel; a crucial phenomenon not considered in other cited 
recent numerical investigations. In fact, in those analy- 
sis, the authors could also have verified the occurrence 
of tetramer states merging to an atom-trimer threshold, 
particularly if they had considered the possibility of a 
four-body independent scaling behavior near a Feshbach 
resonance. This could be done by including, for exam- 
ple, tunable three or four-body potentials, which would 
allow an independent change of the three and four-body 
spectra. The position of the atom-trimer resonance is 
not only a function of the atom-atom scattering and the 
three-body scale, but it also depends on the new four- 
body scale. Our results do not exclude the resonant 
dimer-dimer recombination [M] but also add the possi- 
bility of a resonant atom-trimer recombination. 



Finally, we remark that the four-boson scale can be 
driven near the Feshbach resonance by induced four- 
body forces (coming from the one-channel reduction of 
the atomic interaction) [T5J[57]. Therefore, in this case 
the Efimov ratio percolating the tetramer observables is 
not assured anymore. Other universal scaling functions 
can be derived correlating properties of tetramers, in par- 
ticular the one that correlates the binding energies of 
tetramers attached to different trimers, and as well as by 
extending our framework to the scattering region. The 
possibility that tetramers could be formed and driven in 
cold-atom laboratories, as have been achieved by radio 
frequency association techniques applied successfully to 
measure trimer energies [8], will allow to verify experi- 
mentally the universal scaling relations between tetramer 
properties exhibiting a new scale beyond the trimer one. 
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Appendix A: Momentum space representation of FY 
equations 

In the following we present some details on the FY 
components in the momentum space representation. By 
considering the completeness relations, Eq.(12), between 
the permutation operators, we have 



t\K) = J D 3 u' (u\G tP\u'){u'\K) + J D 3 u' D 3 v (u\G tP\u'){u'\v){v\H) 
+ f D 3 u' D 3 u" (u\G Q tP\u')(u'\P 34 \u")(u"\K), 



(Al) 



ir\H) = I D 3 v' (v\G tP\v')(v'\H) + J D 3 v' D 3 u(v\G tP\v')(v'\u)(u\K) 
D 3 v' D 3 u D 3 v! (v\G tP\v') (v'\u) (u\P 3A \u') (u'\K) . 



r 



Therefore, to evaluate the above coupled equations (Al) By considering that, in the expressions (|A2|) we have to 
we need to obtain the following matrix elements: 



insert the two-body matrix elements (14), and that 



(u\G tP\u') : (v\G tP\v') 



(u\P u \u'), (u'\v), (v'\i 



(A2) 



(A3) 



\P\u') = 



1 6i ui - ni(« 2 ,M2,aO 

dx = '- (A4) 

i K 

5[ u'x - U 1 (u' 2 ,u 2 ,x) ) 5[u 3 -u 3 



II 
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IIi(u2, u' 2 ,X) = 



;U 2 + u 2 



(A5) 



-u| + u' 2 2 + U2U 2 2 X, 



(V\P\V') = *K-"3) *te-*) im 



the matrix elements in (|A2j) are given by the following: 

(u\GotP\u) = Go(wi,M2,lt3) ^ Z) 3 u"(u|t|lt")(u"|P|lt'} 

5 ( u 3 - u 3 



= 47T Go(ui,W2,M3)x(mi)t(£u)- 



(A7) 



^ da;x^ni(u 2 ,«2,a;) 



5 \u'i - Ui{u' 2 ,U2,x) 



and 



where 



(A15) 



5 U 2 - n 6 (« 2 ,«3,x) (5 U 3 - n 7 (v 2 , V3,X 



lie (V2,V 3 ,X) 



n 7 («2,U3,a;) 



V 2 + V 3 



v 2 + v ?i + 2V2V3X, 



V 3 - -V2 



W3 + -1> 2 — U2U3X. 



(A16) 



(A17) 



(v\G tP\v') =G (vx,v 2 ,v 3 ) J D 3 v"{v\t\v") (v"\P\v'), 

= 4nG (vi,V2,v 3 ) x(wi)x(ws)t(£u) (A8) 
5(u 3 - «i) <5(^2 - «4) 



{u\Pzi\u) = 2 / '/.r 



(A9) 



5^u' 2 — n 2 (u 2 ,U3,x)^ <5^u 3 — n 3 (u 2 ,u3,a;)^ 



where 



n 2 (« 2 ,M3,a;) 



II3 («2,W3,a;) 



1 8 
3U2 + -U3 



1 2 64 , 16 

gW 2 + ^ M 3 + ^ u 2«3l, 



U 2 -U3 



,1,2 

U 2 + g"3 - ^2U 3 X 



(A10) 



(All) 



(A12) 



s(v 3 - il 4 (M2,U3,a:)^ (5^2 - IIs(M2,M3.») J 



where 



II4 (M2,M3,x) 



II5 (u2,W3,a:) 



1 2 

2 U2+ 3 U 3 



J«2 + g U 3 + g"2U32; 



u 2 - -u 3 



(A13) 



(A14) 



4 4 



Appendix B: Four-Body Total Wave Functions in 
Momentum Space 

The total four-body wave function, which can be writ- 
ten as 



\ + P + P + Pp)\H), (Bl) 



is composed of eighteen FY components, according to the 
possible arrangements of the four particles (ijkl). We 
have 12 of K— type and 6 of H— type, as follows: 



ijk + I -> 



ijl + k 
ikl + j 
jkl + i 

ij + kl 
ik + jl 
il + jk 



ij ' 




+ l 




K 1 




jk 


f i 


+ 1 




^jk,i 


— Pi 3 Pjk Kij,k 


ki - 




+ 1 




K 1 


= Pik Pjk K ij k 


ij - 


HZ 


+ k 




K k - 


= PkiK i [ k 


Jl- 


-i 


-\-k 




jx- fe 


= Pij Pji i\!j 


li 4 


- j 


+ k 






= P il P j lK i £ l 


ik - 


f / 


+ .] 




K 3 

n ik,l 


= Pjk Pki K ij k 


li 4 


- jfc 


+ .1 




K j 


= PuPikKl l 


kl- 


hi 


+ .J 






= Pik Pik 


jk 




+ i 






= P ij P jk P kl K i l k 


lj- 


- k 


+ i 






= p lk p kj K j i il 


kl- 




+ i 




K ' 


- Pij Pkj '\,Li 


ij 4 


-fc 


+ 1 




Hij,kl 




kl- 


-i 


+ j 




Hkl,ij 


= PikPjl Hij,kl 


ik - 


Hi 


+ 1 




Hik,jl 


= Pjk Hij tk i 




- i ■ 


f k 




Hjl,ik 


= PjkPikPjl Hij,kl 


4 


i ■ 


f k 




H lLjk 


— Pij PjkPikPjl P^ij,kl 


jfc 4- i 


+ 1 




Hjk,u 


ij jk Hij.kl 
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In momentum space, we have 



(«|*) - (u\l\K) + (u\P\K) + (u\P 3 4K) 

+ {u\PP 3A \K) + (u\P 34 P\K) + (u\PP 34 P\K) 
+ (u\l\H) + (u\P\H) + (u\P\H) + (u\PP\H) 
= Y,[Ko + H ]; = l,P,P M ,PP 34 ,...(B2) 

where 

K = {u\0\K) = j D 3 u'(u\0\u')(u'\K), (B3) 
and 

H = (u\0\H) = J D 3 v'{u\0\v'){v'\H) (B4) 

= J D\\u\O\v')G {v[,v' 2 ,v' 3 )x{v' l )U{v' 2 ,v' 3 ). 

In the following we present the matrix elements of Kq 
and Hq explicitly: 

K x = (u\l\K) (B5) 
= K(ui,u 2 ,u 3 ) = G (ui,u 2 ,u 3 ) IC(u 2 ,u 3 ) ; 



(B6) 



K P = (u\P\K) = J dx 12 K\Il P ,n P ,u 3 
= [dx 12 G ( Up,U P ,u 3 J K[U 2 P ,u 3 



where 



np = Up(u 1 ,u 2 ,x 12 ) 



Ui u 2 

2 1 4 2 



19 3 



TT 2 
lip 



^-p(ui,u 2 ,x 12 ) 



Ul - 2 U2 



(B7) 



(B8) 



?i 2 + -u 2 - uiu 2 xi 2 ; 



K P3A = (u\P 3i \K) = X -j dx 23 K(u u Ul 3i ,Up 3 )j (B9) 



dx 23 G [ Ul , Up 3i , u% 4 ) k ( n P34 , n P34 



where 



n 



n 3 
ll p 3 . 



2 _ rr2 



n p 34 ( u 2,U3,a;23) = 



3 U2 + 9 U3 



1 , 64 , 16 

9 U 2 + ^ U 3 + ^fU 2 U 3 X 23 , 



1 



n p 34 (M2, M3,a;23) = 



U 2 - gU 3 



(BIO) 



(Bll) 



-u 2 u 3 x 23 ; 



K P p 3i = (u\PP 34 \K) 



(B12) 



= ^ y dzi2 y rfa;i2, 3 Go^npp 34 ,npp 34 ,npp 34 ^ 



x /c n 2 n 3 



where 



n 



pp 3 . 



n 2 

U PP 3 4 



npp 3 . 



np(wi,-U2,xi 2 ), (B13) 
np 34 (tip^x, u 2 , x 12 ) , u 3 , x 12 ^ , (B14) 

n P 3 4 ( n p( M i, u 2,a;i2),M 3 ,a;i2, 3 ) ; (B15) 



k PmP = (u\p 3i p\K) = iy dx 23 y dz 23il # |n P34 p,np 34 p,np 34 p^ 

= ^ <te23 y ^ 234 Go^np 34 p,n P34 p,np 34 p^/c^np 34 p,np 34 p^, 



(Bi6) 



19 



where 



IIp 34 p = np^Ui,np 34 (u2,U3,X23), 223,1 

np 34 p = np^ui,np 34 (u 2 ,u3, ^23), 223,1^, np 34P = np 34 (u2,u 3 ,22 3 ) ; 



(B17) 



where 



(u\PP 34 P\K) 



dxi2 I dxi23 I ctei2 ; i2,3-K(npp 34 p,npp 34 p,npp 34 p 



= \J rfx i2 J dx i2,3 J ^i2;i2,3G ^npp 34 p,npp 34 p,npp M p^/c^npp 34 p,n 

n PP 3 iP = n P ^ n P ( U l > U 2, X 12 ) , np M ^np («! , U 2 , 2 12 ) , U 3 > ^12,3^ , 2i2;12,3^ , 

n pp 34 p = n P ^ n P ("1^2, 212) , np 34 ^np («i , u 2 , x l2 ) , u 3 , 212,3^ , 2i 2; i 2 ,3^ , 

n pP, 4 P = n p 34 (^ n p( U l ' M 2,2 i2 ), W 3 ,Xi2, 3 



PPmP ' 



(B18) 



(B19) 



(B20) 



(B21) 



For the matrix elements of Hq, we need also to evaluate 
the transformation between the two different representa- 
tions, i.e., the K— and H— type. 



Hi = (u\l\H) = iy da^tf^, (B22) 



where 



where 



^2 v3 



Si = S 2 (u 2 , 1*3,223) 



u 2 + gU 3 



4 4 

M 2 + g W 3 + 3 M 2U 3 223, 



S? = S3(U2,U3,2 23 ) 



1 2 
U 2 U3 

2 3 



4 U 2 + g u 3 - ^u 2 u 3 x 23 ; 



(B23) 



si 



V2 
2jp 



£ 3 
2jp 



where 



nj,(ui,u 2 ,xi2), (B25) 

Si ^np (ui, m 2 , 212) , U3, 212,3^ , (B26) 

Sf^np(ui,w 2 , 212), U3, 212,3^ ; (B27) 

(u\P\H) = ^J d2 23 ff(s p ,S p , Ul ) (B28) 
^ rf2 23 G (s p ,E p , Ul ^S p , Ul 



S p = Sj(w 2 ,U3,22 3 ), 

s p = s 2 (w 2 ,u 3 ,223) ; 



(B29) 
(B30) 



(u\P\H) (B24) 
\ j dx 12 J rf2 12 , 3 H^Sp,Sp,Sp^ 

\ J dx 12 J dx 12 . 3 Go (s P , Sp, Sp^) H (sp, Sp 



pp 



(u|PP|iT) (B31) 

^ y ^212 y d2 12 , 3 ff(^s pp ,s pp ,s pp ^ 
^ y ^12 y d2 12 , 3 G ^s pp ,s pp ,s pp ^ 



x H ( S pp , S pp 
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where 



pp 



Y 2 
^pp 

Y 3 
^pp 



Sf I^IIp (u 1 ,U2,x 12 ) , "3, ^12,3 J , (B32) 

El \U% (u 1 ,u 2 ,xi 2 ) , u 3 , xi2, 3 ^ , (B33) 
U 1 p(u 1 ,U2,x 12 ) . (B34) 



Appendix C: On the numerical approach 

The final set of homogeneous coupled integral equa- 
tions in momentum space, after discretization, defines a 
huge matrix eigenvalue equation, which is solved by it- 
eration within a Lanczos-likc method. This method is 
quite efficient in solving few-body systems. For contin- 
uous momentum and angle variables discretization we 
have used Gaussian-quadrature grid points with hyper- 
bolic and linear mappings, respectively. After consider- 
ing a convenient mesh distribution to have optimal nu- 
merical stability, the number of grid points found neces- 
sary in mapping the momentum variable was set up to 
140. That was achieved by studying the concentration of 
the points in the different relevant regions of momentum 
integration. For more details on the general numerical 
techniques that we are considering, see also Ref. [55] . 

Next, we also address problems concerning to the ac- 
curacy and precision of the results in the unitary limit. 



1. Numerical convergence 

The numerical convergence of our results is exemplified 
with one case, where stability is more difficult to achieve. 
As this happens in the excited states of large scale ratios, 
we choose the case of /i4/ '(13 = 300, in which we can 
verify the existence of up to four tetramers. The number 
of points and, particularly, their distribution, is quite 
critical for the accuracy of the excited-state energies and 
corresponding FY reduced amplitudes. In our systematic 
study to construct the limit cycle shown in Fig. [5j we are 
satisfied to achieve ~1% of deviation from the converged 
results. To be clear, we choose to present in Fig. [13] the 
convergence in respect to the mesh number of points, for 
the first and second excited states. 

As the corresponding FY components become quite 
concentrated at the momentum origin, in this case the 
momentum space discretization Ui is derived from the 
quadrature Gaussian mesh distribution xi = { — 1,+1} 
by the following hyperbolic mapping: 



1 



Ci(l - Xi) + c 2 xC 



(CI) 



with ci = and c 2 — 0.4. In general, we ob- 

serve that for all the cases with smaller scaling ratios, 
or for excited states that are not too close to the contin- 
uum threshold, i.e., about to be unbound, 140 points are 



enough to achieve less than 1% of inaccuracy, within an 
appropriate mesh distribution. 



Iterative Lanczos-type diagonalization algorithm 
for solution of the coupled FY equations 



The coupled FY integral equations (24) can be 



schematically represented as an eigenvalue problem: 

k(E).tp = X(E).ip (C2) 

where the kernel of integral equations k(E) is energy de- 
pendent, \(E) and tp are its eigenvalue and eigenvector, 
correspondingly. The vector i/> is composed of FY com- 
ponents as "0 = {n)- The binding energy of four-body 
bound state can be obtained when one of the eigenval- 
ues obtained from the solution of Eq. (C2| for an input 
energy becomes one (X(E) = 1). To this aim one should 



solve the eigenvalue equation (C2| for a set of input en- 
ergies. 

After discretization of continuous momentum and an- 
gle variables, the kernel of the eigenvalue equation turns 
into a huge matrix with dimensionality of 120 x 120 x 
120 x 40 x 40 ~ 10 9 , where we have used 120 mesh 
points for Jacobi momentum variables (with a hyperbolic 
mapping) and 40 mesh points for angle variables (with 
a linear mapping). Considering that, exact solutions of 
huge matrices by diagonalization are not so efficient, we 
avoid such procedure by using a Lanczos-type technique, 
which is based on iteration. In this way, one can obtain 
the eigenvalues and eigenvectors of a huge matrix, from 
input energies [56J. In the following, we describe some 
details of this technique. 

The iteration procedure is performed with a properly 

starting vector ip Q — ^ , chosen as Gaussian functions 

for both K. and T-L components. After N iterations, one 
obtains the set of vectors {ipi} = ipi,ip 2 , ■■■,ipN where 



= k{E).^_ 



;i = l,2,...,JV 



(C3) 



In the iterative diagonalization approach, an orthonormal 
basis {V'i} can be built up by following a recursive pro- 
cedure and by using the original vector set {ipi}, where 

{i— 1 n 
^-E(^)^f (i = l,2,...,JV). (C4) 
3=1 ' 

In the above, Cj are normalization factors which can 
be obtained by orthogonalization of orthonormal basis 
states ipi-ipj — Sij. By introducing a matrix D, obtained 
from the FY components as 

Dij = ipi.t/jj = (jCiHiY^ J ) 

= KiXj + UiMj = + DS, (C5) 
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where 



Z\j = JCi.tCj 



d,U 2 u\ I du 3 ul lCi(U2,U 3 )lCj(u2,U3), 



= Hi-Uj := J dv 2 v\ J dv 3 vl Hi(v 2 ,v 3 )'Hj(v2,V3), 

(C6) 



the normalization factors Cj can be obtained as: 

1 



(C7) 



The matrix elements ATjj are projection of vectors ipi on 
tpj. i.e. A^. = Vi-V 1 ?) 



Ni^cJDij-J^Nik-NjA 



(C8) 



By expanding the original vector ip, in Eq. (C2|, in 
terms of constructed orthonormal vectors: 



N-l 



(C9) 



and projecting the eigenvalue equation onto tpj, the huge 
dimension of the eigenvalue problem will be reduced to a 
problem of much smaller dimension N — 1 : 



N-l 



(CIO) 



i=l 



where 



By defining and 6^ as coefficients of linear expansion 
of original vectors in terms of orthonormal vectors and 
reversely, 



tpi = ^dijipj, ip i = 'J2b ij i(j j , (C12) 

3=1 3=1 

the matrix elements of M can be obtained as: 

i 

= b *k ak+ij 1 < i, j < N (C13) 



fe=i 



where 




c A-T?k=3+i h ik a kj) J ■ i 



3 



(C14) 
(C15) 

(C16) 



The reduced eigenvalue problem, Eq. (C13l, can be 



solved easily by diagonalization of the matrix M. For a 
given input energy, one can obtain a set of eigenvalues 
and their corresponding eigenvectors. In order to obtain 
the physical binding energy, one should look for an eigen- 
value A = 1 in the obtained eigenvalue spectrum. Clearly, 
the eigenvector corresponding to this eigenvalue gives the 
coefficients gi . Therefore, by using the Eq. ( C9 ) , one can 



(Cll) obtain the physical eigenvector or the FY components. 



[1] K. G. Wilson, Phys. Rev. D 3, 1818 (1971). 

[2] T. Frederico, L. Tomio, A. Delfino, and A. E. A. Amorim, 

Phys. Rev. A 60, R9 (1999). 
[3] P. F. Bedaque, H.-W. Hammer, and U. van Kolck, Phys. 

Rev. Lett. 82, 463 (1999). 
[4] R. F. Mohr, R. J. Furnstahl, H.-W. Hammer, R. J. Perry, 

and K. G. Wilson, Ann. Phys. 321, 225 (2006). 
[5] V. Efimov, Phys. Lett. B 33, 563 (1970); V. Efimov, Sov. 

J. Nucl. Phys. 12, 589 (1971). 
[6] T. Kraemer et al, Nature 440, 315 (2006). 
[7] G. Barontini et al, Phys. Rev. Lett. 103, 043201 (2009). 
[8] T. Lompe et al, Science 330, 940 (2010). 
[9] F. Ferlaino and R. Grimm, Physics 3, 9 (2010). 
[10] R. D. Amado and F. C. Greenwood, Phys. Rev. D 7, 

2517 (1973). 

[11] H. Kroger and R. Perne, Phys. Rev. C 22, 21 (1980). 
[12] S. K. Adhikari and A. C. Fonseca, Phys. Rev. D 24, 416 
(1981). 

[13] H. W. L. Naus and J. A. Tjon, Few-Body Systems 2, 



121 (1987). 

[14] L. Platter, H. W. Hammer, and Ulf-G. Meissner, Phys. 

Rev. A 70, 52101 (2004). 
[15] M. T. Yamashita, L. Tomio, A. Delfino, and T. Frederico, 

Europhys. Lett. 75, 555 (2006). 
[16] H.-W. Hammer and L. Platter, Eur. Phys. J. A 32, 113 

(2007). 

[17] J. von Stecher, J. P. D'Incao, and C. H. Greene, Nature 

Physics 5, 417 (2009). 
[18] Y. Wang and B. D. Esry, Phys. Rev. Lett. 102, 133201 
(2009). 



[19] A. Deltuva, [arXiv:10 09. 1295^1 [physics.atm-clus]; Phys. 

Rev. A 82, 040701(R) (2010). 
[20] J. von Stecher, J. Phys. B: At. Mol. Opt. Phys. 43, 
101002 (2010). 

[21] R. Lazauskas and J. Carbonell, Phys. Rev. A 73, 062717 
(2006). 

[22] F. Ferlaino et al, Phys. Rev. Lett. 102, 140401 (2009). 
[23] M. Zaccanti et al, Nature Phys. 5, 586 (2009). 



22 



[24] S. Pollack et ai, Science 326, 1683 (2009), and 

www.sciencemag.org/cgi / content / full / 1 182840/DC 1 for 

Supporting Online Material. 
[25] G. Modugno, Science 326, 1640 (2009). 
[26] M. R. Hadizadeh, M. T. Yamashita, L. Tomio, A. Delfino, 

and T. Frederico, Phys. Rev. Lett. 107, 135304 (2011). 
[27] S. Nakajima, M. Horikoshi, T. Mukaiyama, P. Naidon, 

and M. Ueda, Phys. Rev. Lett. 105, 023201 (2010); 106, 

143201 (2011). 

[28] L.D. Faddeev, Zh. Eksp. Teor. Fiz. 39, 1459 (1960) [Sov. 

Phys. JETP 12, 1014-1019 (1961)]. 
[29] O.A. Yakubovsky, Yad. Fiz. 5, 1312 (1967) [Sov. J. Nucl. 

Phys. 5, 937 (1967)]. 
[30] S. Weinberg, Phys. Rev. 131, 440 (1963). 
[31] J. A. Tjon, Phys. Lett. B 56, 217 (1975). 
[32] T. Sasakawa, Phys. Rev. C 13, 1801 (1976). 
[33] B. F. Gibson and D. R. Lehman, Phys. Rev. C 14, 685 

(1976); Phys. Rev. C 15, 2257 (1977). 
[34] S. P. Merkuriev, S. L. Yakovlev, and C. Gignoux, Nucl. 

Phys. A 431, 125 (1984); 
[35] N. W. Schellingerhout, J. J. Schut, and L. P. Kok, Phys. 

Rev. C 46, 1192 (1992). 
[36] H. Kamada and W. Glockle, Phys. Lett. B 292, 1 (1992); 

H. Kamada et al., Phys. Rev. C 64, 044001 (2001). 
[37] A. L. Zubarev and V. B. Mandelzweig, Phys. Rev. C 52, 

509 (1995). 

[38] I. N. Filikhin, S. L. Yakovlev, V. A. Roudnev, and B. 

Vlahovic, J. Phys. A: At. Mol. Opt. Phys. 35, 501 (2002). 
[39] R. Lazauskas and J. Carbonell, Few-Body Syst. 34, 105 

(2004). 

[40] S. Bayegan and M. R. Hadizadeh, W. Glockle, Prog. 

Theor. Phys. 120, 887 (2008) 
[41] M. R. Hadizadeh, L. Tomio, and S. Bayegan, Phys. Rev. 



C 83, 054004 (2011). 
[42] R. Perne and W. Sandhas, Phys. Rev. Lett. 39, 788 
(1977). 

[43] S. A. Sofianos, H. Fiedeldey, H. Haberzettl, and W. Sand- 
has, Phys. Rev. C 26, 228 (1982). 

[44] S. Nakaichi-Maeda and T. K. Lim, Phys. Rev. A 28, 692 
(1983). 

[45] A. C. Fonseca, Phys. Rev. C 30, 35 (1984). 
[46] E. O. Alt, P. Grassberger and W. Sandhas, Phys. Rev. 
C 1, 85 (1970). 

[47] T. Frederico, L. Tomio, A. Delfino, M.R. Hadizadeh, 

and M.T. Yamashita, Few-Body Syst. (2011) [DOI 

10.1007/s00601-011-0236-7]. 
[48] L. H. Thomas, Phys. Rev. 47, 903 (1935). 
[49] S. K. Adhikari, A. Delfino, T. Frederico, I. D. Goldman, 

and L. Tomio, Phys. Rev. A 37, 3666 (1988). 
[50] L. Tomio, M.T. Yamashita, T. Frederico, and F. Bringas, 

Laser Physics 21, 1464 (2011). 
[51] M. Th0gersen, D. V. Fedorov, and A. S. Jensen, Phys. 

Rev. A 78, 020501(R) (2008). 
[52] Th. Cornelius and W. Glockle, J. Chem. Phys. 85, 1 

(1986). 

[53] D. Blume and C. Greene, J. Chem. Phys. 112, 8053 
(2000). 

[54] J. P. D'Incao, J. von Stecher, and C. H. Greene, Phys. 

Rev. Lett. 103, 033004 (2009). 
[55] M. R. Hadizadeh and S. Bayegan, Few Body Syst. 40, 

171 (2007). 

[56] A. Stadler, W. Glockle, and P. U. Sauer, Phys. Rev. C 

44, 2319 (1991). 
[57] L. Platter, C. Ji, and D. R. Phillips, Phys. Rev. A 79, 

022702 (2009). 



23 




200 



o 

<n 

E 



3 



150 



100 



50 



1 

(bl) 


AN=2 


1 




\n=i 








N=0 



10 



10 

(in units of \i A 



oa 



o.i 




10 



10 



150 



'5" 









i2 


100- 


'c 




=5 




c 






50- 








0- 




10 



10"" 



10 

(in units of u_ 4 



10- 



CQ. 



0.2- 



0.1 




10 



10 



FIG. 8: (Color on-line) Momentum distribution functions n(ui), n(u2) and n{uz), with scale ratio ^14/^3 = 50, for the ground 
(TV = 0) and the first two (N = 1, 2) excited tetramer levels, normalized as shown in Eq. 30 In the left panels, they are given 



as functions of the corresponding momentum, where all momenta are in units of the momentum scale /J4. In the right panels, 
the momentum distributions are rescaled by the momentum factor /3at 



B4 , where TV is the tetramer level. 
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FIG. 9: (Color on-line) Renormalized momentum distributions as functions of the Jacobi momenta 1^=1,2,3, considering the 
scaling ratios ^14/^3 = 5, when there is only one excited state, and /i4//i3 = 50, when we have two excited states. The units 

are given by the momentum factor Pm = J B^" 1 . For each momentum component Ui, the corresponding distributions of both 
scale ratios, when renormalized to the same maximum value, are shown to reduce almost to the same final form (the labels 
inside the frame are just to indicate the performed calculations). 
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FIG. 10: (Color on-line) The magnitude of the 4B total wave function ^(ui, 112, U3) as a function of Ui and U2, with U3 = 0, for 
the scale ratio [14/(13 = 50. The corresponding contour plots are shown for low momentum region, near the peak of the wave 
function. 




FIG. 11: (Color on-line) As in Fig. 



10 



the wave-function in given in terms of u\ and 113, with 112 = 0. 
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FIG. 13: (Color on-line) Convergence ol numerical results for the first and second excited state energies, in the case that 
= 300. 



